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This  paper  presents  the  assessment  of  stability  domains  for  the  angle  stability 
condition  of  the  power  system  using  Particle  Swarm  Optimization  (PSO) 
technique.  An  efficient  optimization  method  using  PSO  for  synchronizing 
torque  coefficients  and  damping  torque  coefficients  to  identify  the 
angle  stability  condition  on  multi-machine  system.  In  order  to  accelerate  the 
determination  of  angle  stability,  PSO  is  proposed  to  be  implemented  in  this 
study.  The  application  of  the  proposed  algorithm  has  been  justified  as  the 
most  accurate  with  lower  computation  time  as  compared  to  other 
optimization  techniques  such  as  Evolutionary  Programming  (EP)  and 
Artificial  Immune  System  (AIS).  Validation  with  respect  to  eigenvalues 
determination,  Least  Square  (LS)  method  and  minimum  damping  ratio 
confirmed  that  the  proposed  technique  is  feasible  to  solve  the  angle  stability 
problems. 
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1.  INTRODUCTION 

With  the  increase  of  energy  consumption  in  this  age,  a  study  on  the  stability  of  the  power  system 
becomes  a  necessity,  especially  small  signal  stability  analysis  of  power  systems.  This  analysis  is  used  to 
predict  electromagnetic  swing  at  low  frequencies,  as  a  result  of  undisturbed  rotor  swing.  References  [l]-[9] 
have  pointed  out  that  the  stability  of  the  oscillation  in  the  power  system  is  a  very  important  issue.  As  the 
power  system  operation  changes  over  time,  the  stability  of  the  small  signal  in  this  power  system  should  be 
tracked  online.  To  track  the  system,  selected  stability  indicators  are  calculated  from  the  data  provided  over 
time.  These  indicators  are  updated  until  a  constant  value  is  obtained.  In  this  paper,  synchronizing  torque 
coefficient,  A,  and  damping  torque  coefficient,  are  used  as  stability  indicators.  Eor  a  system  to  achieve  a 
stable  condition,  both  A,  and  values  must  be  positive  [5] -[8]. 

The  Least  Square  method  (LS)  is  one  of  the  techniques  in  finding  this  A,  and  value,  which  has 
been  used  as  a  static  parameter  estimator  [6] -[8].  However,  long  calculation  times  and  data  updating 
requirements  are  the  weakness  in  the  LS  method.  This  method  also  requires  monitoring  throughout  the 
duration  of  the  swing. 

Computational  intelligence  techniques  have  been  widely  used  in  solving  power  system  stability 
problems.  One  of  them  is  Evolutionary  Programming  (EP).  This  algorithm  is  heuristic  population-based 
search  methods  that  used  both  random  variation  and  selection.  In  this  paper,  new  technique  called  PSO  is 
proposed.  It  brought  about  the  performance  beyond  EP  method  in  searching  the  optimal  solution  with  faster 
computation  time. 
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This  paper  proposes  an  efficient  technique  for  estimating  synchronizing  and  damping  torque 
coefficients  in  solving  angle  stability  problems.  This  technique  is  used  to  estimate  the  value  of  and 
from  the  machine  responses  namely,  the  change  in  rotor  angle,  Ad(t),  the  change  in  rotor  speed,  A(o(t)  and  the 
change  in  electromechanical  torque,  ATJt).  The  goal  is  to  minimize  the  error  of  the  estimated  coefficients. 
The  online  estimation  technique  for  synchronizing  and  damping  torque  coefficients  has  been  tested  in  IEEE 
9-Bus  system. 

This  study  presents  PSO  based  technique  for  angle  stability  assessment  in  multi-machine  system. 
Through  this  study,  a  mathematical  model  for  multi-machine  for  the  assessment  of  angle  stability  is 
developed.  PSO  is  used  to  optimize  the  objective  function,  J  with  and  Kd  as  the  control  variables.  Once  the 
J  value  has  been  maximized,  and  are  analyzed  which  determines  whether  the  rotor  angle  is  stable  or 
vice  versa.  The  performance  of  PSO  is  then  compared  to  EP  and  AIS .  Results  obtained  from  the  experiment 
were  then  verified  with  minimum  damping  ratio,  and  eigenvalues,  1. 


2.  PROBLEM  FORMULATION 

To  demonstrate  the  potential  of  proposed  technique  in  angle  stability  assessment  for  multi-machine 
system,  IEEE  9-Bus  System  model  has  been  selected.  Three  generators  called  Generator  1,  Generator  2  and 
Generator  3  are  connected  to  the  buses  named  Bus  1,  Bus  2  and  Bus  3,  respectively.  The  single  line  diagram 
of  the  test  system  is  showed  in  Eigure  1.  The  parameters  of  the  system  are  shown  in  [9]. 


6  7  8 


Eigure  1.  Single  line  diagram  for  IEEE  9-bus  system 


2.1.  Proposed  Philips-Heffron  Model  for  Multi-machine  System 

A  proposed  Phillips -Heffron  model  for  the  multi  machine  system  is  developed  and  showed  in  Eigure 
2.  It  is  developed  based  on  the  single  machine  of  Philips-Heffron  model  [10]. 
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Eigure  2.  Phillips-Heffron  model  for  multi-machine  system 


Kd  is  the  damping  torque  coefficient,  H  is  the  inertia  constant,  and  Ta  are  the  circuit  constant  and 
time  constant  of  the  exciter  oscillation  system,  respectively,  coo  is  equal  to  lirfo-  Ks  is  a  function  of  the  ratio  of 
impedance.  Kj,  K2,  T3,  K4,  Ks  and  are  constants  which  consist  of  function  related  to  the  operating  real  and 
reactive  loading,  electrical  torque,  rotor  speed,  rotor  angle  as  well  as  the  excitation  levels  in  the  generator. 

2.2.  The  Developed  Mathematical  Modelling 

The  mathematical  modelling  can  be  derived  for  the  proposed  Phillips-Heffron  model  for  multi 
machine  system  of  Eigure  2.  They  are  presented  in  the  following  mathematical  equations: 
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Details  on  Equation  1  to  5  are  shown  in  [10].  The  Equation  1  to  5  can  be  rewrite  into  matrix  form  as 
followed: 
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Xi  and  Ui  are  the  state  vector  and  input  signal  vectors  for  i  generators,  respectively.  The  system 
matrix  A,  is  a  function  of  the  system  parameters  with  i  generators  that  depends  on  the  opening  conditions. 
The  perturbation  matrix  Bi  depends  on  the  system  parameters  only. 

2.3.  Synchronizing  and  Damping  Torque  Coefficients 

The  correlation  between  the  change  in  estimated  electromagnetic  torque  deviation,  AT^Jt)  with  the 
change  in  rotor  angle,  Ad/f]  and  the  change  in  rotor  speed,  AcOi(t)  for  the  generator  can  be  expressed  as: 

Ar,,,  (f)  =  (0  +  =  (9) 

where  fCi  and  are  synchronizing  torque  coefficient  and  damping  torque  coefficient  for  the  generator, 
respectively,  m  is  the  number  of  generators. 

The  justification  of  the  stability  of  a  linear  system  can  be  done  by  the  estimation  of  and  K^. 
Positive  value  of  both  and  will  validate  the  system  as  stable.  If  the  system  has  positive  and  negative 
K^,  the  system  is  defined  as  in  oscillatory  instability  condition.  On  the  other  hand,  if  and  respectively 
show  negative  and  positive  value,  the  system  is  considered  as  non-oscillatory  instability  condition.  In 
general,  the  system  is  said  to  be  unstable  if  either  one  of  the  torque  coefficients  is  negative. 

The  stability  evaluation  of  a  linear  system  can  be  predicted  by  referring  to  and  values.  A  stable 
system  is  guaranteed  if  both  and  values  are  positive.  Eigure  3  illustrated  a  stable  angle  stability  resulted 
from  a  positive  value  of  both  and  fC;  [  1 1  ■ 


A<5 


t 


Eigure  3.  Complex  plane  of  Ar/tj  and  A8(t)  response  for  stable  condition 


If  the  linear  system  has  positive  and  negative  K^,  the  system  is  defined  as  in  oscillatory  instability 
condition,  which  is  due  to  lack  of  adequate  damping  torque.  The  effect  of  the  oscillatory  instability  condition 
can  be  detected  from  the  increment  of  amplitude  oscillations  of  the  rotor.  Eigure  4  illustrated  unstable 
conditions  for  angle  stability  resulted  from  positive  value  of  and  negative  value  of  K^. 
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Eigure  4.  Complex  plane  of  ATJt)  and  Ad(t)  response  for  oscillatory  instability  condition 


Non-oscillatory  instability  occurred  if  and  K^,  respectively  show  negative  and  positive  value.  This 
is  due  to  the  absence  of  automatic  voltage  regulators,  resulting  lack  of  sufficient  synchronizing  torque.  This 
condition  can  be  verified  from  steady  increment  of  rotor  angle  response.  Eigure  5  illustrated  unstable 
conditions  for  angle  stability  resulted  from  negative  value  of  and  positive  value  of  Kj. 


Eigure  5.  Complex  plane  of  ATJt)  and  Ad(t)  response  for  non-oscillatory  instability  condition 


2.4.  Eigenvalues,  Damping  Ratio  and  Damping  Factor 

The  scalar  parameter  of  eigenvalues,  1  can  be  derived  as  follows  [1]: 

(A-/IIV  =  0  (10) 

where  A  is  an  n  x  n  matrix  and  ^  is  an  n  x  1  vector. 

The  n  solutions  of  1=17,  ^2,  ■■■,  K  are  eigenvalues  of  A.  The  eigenvalue  can  be  stated  as  follows: 


T;  =  0-,.  ±  Jftt;  (11) 

where  cr,  is  the  real  part  of  the  eigenvalue  and  co,  is  the  imaginary  part  of  the  eigenvalue. 

The  linear  system  is  stable  if  all  eigenvalues  have  negative  real  parts.  The  damping  ratio,  for  the 

eigenvalue  is  defined  as: 

=  -<^il  (12) 

The  linear  system  is  certainly  in  stable  condition  if  all  damping  ratio  have  positive  value.  Eor 
simplification  purposes,  only  the  minimum  value  of  damping  ratio,  for  the  linear  system  is  selected  to 
verify  the  result. 

2.5.  Least  Square  (LS)  Method 

Least  square  (LS)  technique  is  used  to  minimize  the  sum  of  the  square  of  the  differences  between 
ATe(t)  and  ATe,(t).  The  error  is  defined  as  [5],  [7]: 

E{t)=ATXt)-ATjt)  (13) 


where  ATJt):  the  electrical  torque  AEJt):  the  estimated  electrical  torque. 
The  value  of  t,oiai  is  as  follows: 
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t,otai=NxT  (14) 

where  N  is  the  number  of  samples  T  is  the  sampling  period. 

To  estimate  the  correct  value  for  K,  and  K^,  the  period  of  t,o,ai  must  be  chosen.  Different  value  of 
will  result  an  inaccurate  value  for  and  K^.  [7]  and  [8]  have  stated  that,  the  suitable  value  for  t,otai  that 
makes  and  constant  during  the  oscillation  period  is  the  value  of  the  entire  period  of  oscillation.  In 
matrix  notation,  the  above  problem  can  be  described  by  an  over -determined  system  of  linear  equations  as 
follows: 


^TXt)  =  ^T^Xt)+E{t)=Cx  +  E{t)  (15) 

C  =  [Aj(f)  A®(t)]  (16) 

x  =  [k,  kJ  (17) 

where  /ITJt):  the  electrical  torquezir^/t):  the  estimated  electrical  torqueiiff):  differences  (error)  between 
ATJt)  and/17),/f) 

The  estimated  vector,  x  is  such  that  the  function,  F(x)  is  minimized,  where 

T’(x)  =  [AT^  -  AxJ  ■  [at;  -  Ax]  (1 8) 

In  this  case,  x  will  be  given  by: 

+  X  =  {a'^  ■  aY  ■  A'^  ■  AT^  =  A‘ ■  AT^  (19) 

where  A‘  is  the  left  pseudo  inverse  matrix. 


By  solving  Equation  19,  it  gives  the  values  of  and  for  the  corresponding  operating  point. 
Although  the  calculated  values  are  accurate,  the  application  of  LS  method  is  time  consuming  and  requires  the 
entire  period  of  oscillation  [7] -[8].  Due  to  that,  new  indicator  is  needed. 


3.  COMPUTATIONAL  INTELLIGENCE  APPROACHES 

Lately,  the  use  of  Artificial  Intelligence  (AI)  technology  is  synonymous  in  solving  power  system 
problems.  AI  techniques  utilize  the  logic  and  knowledge  representations  of  expert  systems,  artificial  neural 
network  (ANN)  [11]-[12]  and  evolutionary  computation  (EC).  The  EC  field  includes  PSO  [13]-[17], 
Evolutionary  Programming  (EP)  [8],  [18]-[19],  Artificial  Immune  Systems  (AIS)  [20]-[22]  and  Eirefly 
Algorithm  (PA)  [23]-[24].  In  this  study,  EP,  PSO  and  AIS  were  selected  as  optimization  techniques. 

3.1.  Particle  Swarm  Optimization 

PSO  was  introduced  by  Dr.  Russ  Eberhart  and  Dr.  James  Kennedy  in  1995.  PSO  is  an  evolutionary 
based  optimization  technique,  which  imitates  the  behaviour  of  birds  flocking  and  fish  schooling.  The  PSO 
algorithm  is  started  with  initialization,  followed  by  the  update  of  velocity  and  position,  fitness  calculation,  the 
best  position  update  and  convergence  test.  The  flow  chart  which  represents  the  PSO  algorithm  is  illustrated  in 
Pigure  6. 
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Eigure  6.  Elowchart  for  PSO  algorithm 


Parameters  such  as  acceleration  coefficients  and  inertia  weight  determine  the  performance  of  PSO 
technique.  The  inertia  weight  carries  a  huge  impact  on  the  convergence  capability  of  PSO.  A  small  inertia 
weight  can  facilitate  a  local  exploration,  while  a  large  inertia  weight  will  facilitate  an  exploration  globally.  In 
this  study,  the  evaluation  of  inertia  weight  is  performed  using  Linearly  Decreasing  Inertia  Weight  (LDIW) 
method.  In  the  LDIW  approach,  the  inertia  weight,  co  is  designed  to  decrease  linearly  during  the  evolution  of 
the  PSO,  from  a  maximum  value,  (Omax,  decreases  to  a  minimum  value,  cOniin  [17].  The  equation  of  inertia 
weight  using  LDIW  approach  is  given  in: 

(Ok  =  <«max  -  {k/ )  (20) 

where  iter^^  is  the  maximum  number  of  iteration  and  k  is  the  current  iteration. 

3.2.  Evolutionary  Programming 

The  Evolutionary  Programming  (EP)  uses  the  models  of  biological  evolutionary  process  to  obtain 
the  solution  for  complex  engineering  problems.  Invented  by  D.  Eogel  in  1962,  the  optimal  solution  search 
method  using  EP  technique  is  accomplished  in  a  parallel  method  within  the  parameter  search  space.  The 
process  of  Evolutionary  Programming  (EP)  algorithm  includes  initialization,  statistical  evaluation,  fitness 
calculation,  mutation,  combination  and  selection.  The  overall  process  of  EP  algorithm  is  given  in  [8]. 

3.3.  Artificial  Immune  Systems 

Artificial  Immune  System  (AIS)  and  EP  share  many  common  aspects  of  optimization  techniques. 
EP  is  based  on  the  natural  evolution  model,  while  AIS  tries  to  benefit  from  the  characteristics  of  a  human 
immune  system.  AIS  involves  initialization,  statistical  evaluation,  fitness  calculation,  cloning,  mutation, 
combination  and  selection.  The  whole  process  is  given  in  [20]. 

3.4.  Objective  Functions 

In  this  study,  the  objective  function  formulated  is  based  on  the  differences  of  the  electromagnetic 
torque  and  the  estimated  electromagnetic  torque  of  the  generator,  ATJt)  and  AT^Jt),  respectively,  as 
shown  in  Equation  21.  This  difference  or  error  is  estimated  for  calculating  A,  and  for  every  generator  in 
the  system.  PSO  optimization  technique  is  used  to  minimize  the  error  with  and  being  the  control 
variables  [7]. 
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J  i  =  inv\ 


1-F 


Ar,,(0-Ar,,„.(0 


(21) 


where  m  is  the  number  of  generators.  Hence,  the  objective  function  can  be  defined  as: 

Maximize  (7,) 

From  the  optimized  J  value,  decision  can  be  made  to  identify  the  angle  stability  based  on  the  values 
of  Ks  and  K^. 

3.5.  Algorithm  for  Angle  Stability  Assessment 

To  estimate  successfully  the  maximum  value  of  7„  the  calculation  process  of  Aj,  and  for  the 
generator  is  conducted  repeatedly.  The  following  process  is  implemented: 

a.  Calculate using /Id/t),  Jco/t)  and  the  estimated  torque  coefficients  using  Equation  9. 

b.  Evaluate  7,  using  Equation  21 . 

c.  If  7,  is  smaller  than  1.00,  vary  the  value  of  K,/  and  Kj/  and  repeat  step  (i)  and  (ii)  with  new  generated 
^S/t)  and  AcOi(t)  sample  data  until  7,  reaches  1.00  or  all  sample  data  were  used. 

The  process  is  simplified  in  the  form  of  a  flowchart,  as  shown  in  Eigure  7. 


Eigure  7.  Elowchart  for  and  Kj  estimation  process 


Table  1  tabulates  the  parameters  used  in  PSO  optimization  process.  There  are  four  parameters 
needed  for  the  PSO  as  the  optimization  technique  i.e.  ci  and  C2  are  acceleration  coefficients,  co„,„  and  cOmax  are 
the  minimum  inertia  weight  and  the  maximum  inertia  weight,  respectively. 


Table  1.  Parameters  of  EP,  AIS  and  PSO 


Parameters 

EP 

AIS 

PSO 

C]  and  C2 

- 

- 

0.9 

COmin 

- 

- 

0.04 

COmax 

- 

- 

0.09 

P 

0.05 

0.05 

- 

Population  Size 

20 

20 

20 
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4.  RESULTS  AND  DISCUSSION 

In  this  study,  the  performance  of  PSO  technique  in  estimating  Ks  and  are  conducted  on  IEEE  9- 
bus  system.  Three  samples  of  data  of  rotor  angle,  ASft),  rotor  speed,  Aoj(t)  and  electrical  torque,  AT^(t)  for  all 
three  generators  are  produced  in  Matlab  Simulink  environment.  To  simulate  various  stability  cases,  three 
different  value  of  reactive  load  at  Bus  5  is  used.  The  values  of  the  reactive  load  at  Bus  5  are  chosen  in  such  a 
way  that  three  scenarios  can  be  emulated,  namely  unstable,  critically  stable  and  stable  conditions  as  tabulated 
in  Table  2. 


Table  2.  Three  Different  Loading  Conditions 


Cases 

Reactive  Load  at  Bus  5 

Case  1  (unstable  condition) 

160  MVar 

Case  2  (critically  stable  condition) 

60  MVar 

Case  3  (stable  condition) 

10  MVar 

Responses  of  Ad(t),  A(o(t)  and  ATJt)  for  Cases  1,  2  and  3  are  shown  in  Eigure  8,  Eigure  9  and  Eigure 
10,  respectively.  To  evaluate  the  capability  of  the  proposed  scheme,  the  minimum  damping  ratio  and 
eigenvalues  1  are  calculated  and  were  set  as  a  benchmark.  In  stable  condition,  S„,i„  will  give  a  positive  result 
meanwhile  all  values  of  1  will  show  negative  values.  On  the  contrary,  if  the  system  is  unstable,  it  will  result 
to  negative  whereas  at  least  one  of  1  will  give  positive  value.  Erom  the  result  of  1  and  first  case  is 
unstable  and  the  last  two  cases  are  stable.  Eigure  8  illustrates  the  responses  for  unstable  condition  for  all 
generators  in  Case  1 . 

In  Eigure  8(a),  the  response  shows  the  angle  deviation  AS(t)  of  Generator  1,  Generator  2  and 
Generator  3  for  Case  1.  Erom  the  response,  it  clearly  shows  that  the  response  damping  for  all  generators  is 
increasing  dramatically  since  the  simulation  starts.  The  same  phenomenon  can  be  seen  in  speed  deviation 
Aco(t)  and  torque  deviation  ATJt)  for  Case  1  as  shown  in  Eigure  8(b)  and  Eigure  8(c),  respectively.  The 
damping  increment  of  the  response  for  all  three  figures  indicates  that  Case  1  is  an  unstable  case. 


(a)  Angle  deviation 


(b)  Speed  deviation 


(c)  Torque  deviation 

Eigure  8.  Responses  for  unstable  condition  for  all  generators  in  Case  1 
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Table  3  tabulates  the  comparison  of  K„  Kj,  fitness  J  and  computation  time  for  three  different 
generators  from  EP,  AIS,  PSO  and  LS  method  for  Case  1.  In  this  study,  no  fitness  value  calculated  for  LS 
method  as  it  is  not  an  optimization  technique.  From  Table  III,  PSO,  EP,  AIS  and  LS  give  negative  values  of 
Ks  and  for  first  generator  Gl.  For  second  and  third  generators,  G2  and  G3,  all  methods  give  positive 
and  negative  K^.  In  this  paper,  value  calculated  with  LS  method  was  selected  as  a  benchmark  and  all  value 
calculated  using  EP,  AIS  and  PSO  estimation  techniques  were  compared  with  this  value.  In  Table  III,  the 
estimated  value  of  and  using  EP,  AIS  and  PSO  showed  similar  values  with  the  values  estimated  by  LS. 
Among  them,  PSO  gives  the  closest  value  compared  to  LS  for  all  generators,  followed  by  EP  and  finally  AIS 
method.  It  shows  that  PSO  technique  manages  to  bring  the  most  accurate  value  of  and  compared  to  EP 
and  AIS.  Since  the  values  of  and  for  all  three  generators  are  negative,  this  result  indicates  that  all 
generators  in  Case  1  are  unstable. 


Table  3.  Comparisons  of  EP,  AIS,  PSO  and  LS  Method  for  Case  1 


Gen. 

Tech 

Ks 

K, 

J 

Computation  Time 

EP 

-0.2563 

-3.2631 

0.8901 

18.0  s 

1 

AIS 

-0.2451 

-4.9650 

0.8561 

7.5  s 

PSO 

-0.2769 

-2.4385 

1. 0000 

20.8  s 

LS 

-0.2638 

-2.9425 

- 

- 

EP 

2.3026 

-4.0787 

0.8106 

19.2  s 

A 

AIS 

2.5156 

-3.3527 

0.8071 

7.5  s 

Z 

PSO 

2.I7I7 

-5.7955 

1. 0000 

19.5  s 

LS 

2.2533 

-5.0252 

- 

- 

EP 

I.4I69 

-1.963 

0.8906 

19.2  s 

O 

AIS 

1. 3 148 

-I.9I7I 

0.8891 

7.5  s 

J 

PSO 

1.4909 

-2.0073 

1. 0000 

19.5  s 

LS 

1.4572 

-1.9558 

- 

- 

From  all  methods,  only  PSO  manage  to  achieve  fitness  value  of  1.000  for  all  generators.  The  second 
highest  of  fitness  value  is  EP,  with  range  between  0.8506  ~  0.8906.  AIS  give  the  lowest  value  of  fitness,  with 
range  between  0.8071  ~  0.8891.  In  terms  of  computation  time,  AIS  is  the  fastest  technique,  manage  to  finish 
the  optimization  process  in  7.5  seconds.  PSO  and  EP  are  in  the  same  level,  able  to  finish  the  optimization 
process  between  18-21  seconds.  LS  technique  is  not  an  optimization  method,  so  no  iteration  value  is 
recorded. 

Table  4  shows  the  results  of  X  and  for  Case  1.  From  eigenvalues  X  point  of  view,  one  of  the 
eigenvalues  is  0.0512  which  is  positive  value.  This  confirms  that  Case  1  is  an  unstable  case.  Negative  value 
of  minimum  damping  ratio  was  also  verified  the  same  conclusion. 


Table  4.  The  results  of  Eigenvalues,  Minimum  Damping  Ratio  for  Case  1 


(min 

-1.0 

0.1098+jl4.7919,  0.2803+j7.4378,  0.0512,  -0.0003, 
-31.866,  -26.8641+j7.472,  -24.1887+j  2.653,  -17.892. 

Figure  9(a),  Figure  9(b)  and  Figure  9(c)  show  the  responses  of  A8(t),  Aa)(t)  and  ATJt)  for  Case  2, 
respectively.  In  Figure  9(a),  the  damping  of  angle  deviation  for  all  generators  is  decreasing  gradually.  When 
the  simulation  time  reaches  50  s,  damping  for  Gl,  G2  and  G3  are  minimum  but  still  not  fully  damped. 
Among  three  generators,  G2  and  G3  show  almost  the  same  responses.  On  the  other  hand,  response  of  the  first 
generator  Gl  is  the  most  fluctuate.  The  same  response  patent  can  be  seen  in  Figure  9(b)  and  Figure  9(c), 
which  show  speed  deviation  and  torque  deviation  of  Case  2,  respectively.  Same  as  Figure  9(a),  the  response 
of  Gl  is  the  most  fluctuate  compared  to  G2  and  G3  for  both  Figure  9(b)  and  Figure  9(c).  These  results 
indicate  that  Case  2  is  a  stable  case. 
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(a)  Angle  deviation 


(b)  Speed  deviation 


(c)  Torque  deviation 

Figure  9.  Responses  for  stable  condition  for  all  generators  in  Case  2 


The  value  of  torque  coefficients  K,  and  together  with  result  of  fitness  J  and  computation  time  for 
Case  2  are  tabulated  in  Table  5.  For  Case  2,  EP,  AIS  and  PSO  method  give  identical  positive  values  of  both 
Ks  and  K^,  which  verify  that  it  is  a  stable  case.  Despite  it  is  a  stable  case,  LS  method  has  demonstrated 
different  results,  a  negative  value  of  for  generator  G1  and  negative  value  for  both  and  for  G3.  This 
shows  that  LS  technique  failed  to  provide  accurate  results  involving  difficult  cases  such  as  Case  2.  The 
results  of  fitness  J  and  computation  time  shown  in  Table  V  are  almost  consistent  with  the  previous  Case  1. 
PSO  score  1.000  in  fitness  value  for  all  generators.  AIS  is  the  fastest  computation  technique,  which  acquired 
4. 5-7. 5  seconds  to  finish  the  optimize  process.  EP  become  the  worst  of  all  three  method,  with  the  slowest 
simulation  takes  24  seconds  to  converge.  Overall,  PSO  is  the  best  technique  to  achieve  the  accurate 
assessment  in  accepted  iteration  limit. 


Table  5.  Comparisons  of  EP,  AIS,  PSO  and  LS  method  for  Case  2 


Gen. 

Tech. 

Ks 

K, 

J 

Computation  Time 

EP 

2.5026 

3.0787 

0.8106 

19.2  s 

1 

AIS 

2.5156 

3.3527 

0.8071 

7.5  s 

PSO 

2.1717 

2.7955 

1.0000 

19.5  s 

LS 

1.2090 

-0.7011 

- 

- 

EP 

2.1959 

6.9951 

0.8695 

24.0  s 

AIS 

2.2500 

7.4038 

0.8420 

4.5  s 

PSO 

1.8114 

6.5481 

1.0000 

19.5  s 

LS 

0.9002 

4.5534 

- 

- 

EP 

0.9973 

0.0316 

0.7870 

20.4  s 

AIS 

0.9703 

0.0103 

0.7841 

7.5  s 

J 

PSO 

1.8537 

0.0340 

1.0000 

18.2  s 

LS 

-0.0012 

-0.0312 

- 

- 
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Table  6  tabulates  the  value  and  1  for  Case  2.  The  value  of  Cmm  is  0.0065  which  is  nearly  closed 
to  negative  value.  For  eigenvalues,  although  all  values  are  negative,  one  of  the  1  values  is  almost  zero  that  is 
near  to  positive  value.  This  shows  that  Case  2  is  considered  as  stable  case  even  though  it  is  almost  near  to 
unstable  region. 


Table  6.  The  Results  of  Eigenvalues,  Minimum  Damping  Ratio  for  Case  2 


0.0065 

-0.3103  +jl4.348,  -0.0730  +jl  1.282,  -0.0012,  -0.0001, 
-32.591,  -30.772,  -25.589  +j8.514,  -18.474,  -17.2905. 

The  responses  of  Aco(t)  and  ATJt)  for  Case  3  are  shown  respectively  in  Figure  10(a),  Figure 
10(b)  and  Figure  10(c).  In  Figure  10(a),  the  angle  deviation  for  G1  is  completely  damped  about  15  s  after  the 
simulation  started.  It  is  follows  by  Ad(t)  of  G3  which  stop  fluctuate  about  18  s,  and  G2  at  the  23^''  seconds 
after  the  simulation  started.  The  speed  deviation  and  torque  deviation  for  Case  3  also  show  the  same  response 
patent,  as  shown  in  Figure  10(b)  and  Figure  10(c),  respectively.  Based  on  these  three  types  of  responses. 
Case  3  is  classified  as  one  of  stable  cases.  More  than  that.  Case  3  is  damp  faster  than  Case  2. 


(a)  Angle  deviation  (b)  Speed  deviation 


(c)  Torque  deviation 

Figure  10.  Responses  for  stable  condition  for  all  generators  in  Case  3 


Table  7  shows  the  comparisons  of  K^,  Kj,  fitness  J  and  computation  time  for  Case  3.  PSO,  EP,  AIS 
and  FS  estimation  techniques  give  positive  values  of  and  Kj  for  all  generators  Gl,  G2  and  G3.  This  result 
indicates  that  Case  3  is  stable  case.  In  this  case,  PSO  demonstrated  the  most  identical  value  of  and  to 
FS  method  compared  to  the  other  two  techniques.  Also,  from  the  result,  PSO  can  calculate  the  fitness  with 
1.000  scores  for  all  three  generators.  On  the  other  hand,  EP  and  AIS  calculate  the  highest  fitness  for  G2  with 
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0.8534  and  0.8067  scores,  respectively.  Although  AIS  is  the  fastest  optimization  technique,  the  computation 
time  for  EP  and  PSO  are  still  acceptable. 


Table  7.  Comparisons  of  EP,  AIS,  PSO  and  LS  method  for  Case  3 


Gen. 

Tech. 

Ks 

K.I 

J 

Computation  Time 

EP 

6.7984 

3.6458 

0.8157 

19.2  s 

1 

AIS 

7.5737 

3.8709 

0.8051 

4.5  s 

1 

PSO 

5.1308 

3.4876 

1.0000 

19.5  s 

LS 

4.0870 

4.2304 

- 

- 

EP 

5.1562 

11.1108 

0.8534 

18.0  s 

AIS 

6.1249 

14.9902 

0.8067 

7.5  s 

L 

PSO 

4.4353 

10.3909 

1.0000 

20.8  s 

LS 

2.2301 

7.8529 

- 

- 

EP 

2.6436 

8.1310 

0.7939 

13.2  s 

"5 

AIS 

2.6122 

8.0070 

0.7871 

7.5  s 

PSO 

2.8125 

9.0231 

1.0000 

15.6  s 

LS 

1.0323 

7.1655 

- 

- 

The  result  of  1  and  for  Case  3  are  showed  in  Table  8.  Positive  value  of  and  all  negative 
value  of  eigenvalues,  1  verify  the  same  conclusion. 


Table  8.  The  results  of  Eigenvalues,  Minimum  Damping  Ratio  for  Case  3 


1 

0.0499 

-0.2026+j  14.3629, -0.3413+j  10.8919, 
-0.000 l+j0.0020,  -32.5934,  -30.3290, 
-25.5359+j  8.1913,  -18.9565,  -17.0365. 

5.  CONCLUSIONS 

This  study  has  presented  the  capability  of  the  proposed  PSO  technique  in  the  assessment  of  angle 
stability,  specifically  in  multi  machine  system.  IEEE  9-bus  test  system  is  chosen  for  the  study.  PSO  method  is 
proven  as  the  most  reliable  optimization  technique  compared  to  EP  and  AIS  method,  although  all  three 
techniques  managed  to  calculate  correctly  all  cases  whether  it  is  stable  or  unstable.  Optimization  by  PSO  has 
achieved  the  highest  accuracy  compared  to  the  other  two  techniques.  PSO  produced  the  highest  accuracy  for 
all  cases,  compared  to  EP  and  AIS.  Prom  the  iteration  perspective,  AIS  recorded  the  fastest  simulation  time, 
while  PSO  and  EP  were  almost  the  same.  Despite  this  disadvantage,  the  time  consumed  for  PSO  simulation 
process  is  still  minimal  and  acceptable. 
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